Quantum chaos in optical systems: The annular billiard 
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We study the dielectric annular billiard as a quantum chaotic model of a micro-optical resonator. 
It differs from conventional billiards with hard- wall boundary conditions in that it is partially open 
and composed of two dielectric media with different refractive indices. The interplay of reflection 
and transmission at the different interfaces gives rise to rich dynamics of classical light rays and 
to a variety of wave phenomena. We study the ray propagation in terms of Poincare surfaces of 
section and complement it with full numerical solutions of the corresponding wave equations. We 
introduce and develop an S-matrix approach to open optical cavities which proves very suitable 
for the identification of resonances of intermediate width that will be most important in future 
applications like optical communication devices. We show that the Husimi representation is a useful 
tool in characterizing resonances and establish the ray-wave correspondence in real and phase space. 
While the simple ray picture provides a good qualitative description of certain system classes, only 
the wave description reveals the quantitative details. 



PACS numbers: 05.45-Mt, 03.65.Sq, 42.25.-p, 42.60.Da 



I. INTRODUCTION 

Billiard systems of many kinds have proven to be very 
fruitful model systems in the field of quantum chaos. 
The methods of investigation are well established both 
for the classical dynamics and for the quantum mechani- 
cal behaviour, with semiclassical methods describing the 
transition from quantum to classical properties. With 
the growing interest in quantum chaos and in meso- 
scopic physics, new systems have entered the stage, in- 
cluding systems exhibiting chaos of classical waves such 
as (macroscopic) microwave billiards [1,2], acoustic res- 
onators [3] as well as deformed microcavities [5-10] which 
can operate as microlasers [11,12]. To describe these 
(two-dimensional) systems one can exploit the analogy 
between the stationary Schrodinger equation and the 
Helmholtz equation for (classical) waves [13]. Quantum 
chaotic experiments using microwave cavities or other 
classical waves (e.g., acoustic or water waves) are based 
on this mathematical equivalence, sec [14] for a review. 
Most of the investigated systems are hard- wall billiards. 
However, for the class of optical, or dielectric, model sys- 
tems the billiard boundary manifests itself by a change in 
the index of refraction allowing for reflection and trans- 
mission of light. The limit of closed systems is ap- 
proached as the difference in the refractive indices reaches 
infinity. 

We emphasize that the openness of optical systems ex- 
tents the set of interesting questions with respect to those 
for closed billiards. In this paper we suggest a further 
extension of the class of open optical cavities by consid- 
ering two regions with different refractive indices inside 
the cavity, which leads to an additional refractive inter- 
face between the two dielectrics inside the resonator. The 
interplay between refraction inside the billiard and par- 
tial reflection at the outer billiard boundary gives rise to 



a variety of phenomena in the classical ray dynamics and 
correspondingly in the wave description of such systems. 

The model we study is the annular billiard shown in 
Fig. 1. It consists of a small disk of radius i?2 placed 
inside a larger disk of radius R\ with a displacement 
5 of the disk centres. This system is well-known from 
quantum and wave mechanical studies of the hard-wall 
configuration [15-17] with non-vanishing wave-functions 
only in the annular region. It carries features of a ray- 
splitting system [18] when each disk is characterized by a 
(stepwise constant) potential (unlike the situation we will 
discuss here, see Sec. II B). Here we consider disks char- 
acterized by indices of refraction, n\ and ri2, respectively, 
with the index of the environment fixed at no = 1 [19]. 
We will study billiard materials with rii > 1 such that 
confinement by total internal reflection is possible. Then 
methods well-known from the description of classical dy- 
namical systems, such as the use of Poincare's surfaces 
of section, can be employed to describe the ray dynam- 
ics. Note that whispering gallery modes in the dielectric 
annular billiard with a metallic inner disk have been dis- 
cussed in Ref. [20]. A detailed study of periodic orbits 
in a specific hard-wall configuration, together with the 
expected consequences on the electromagnetic scattering 
problem was performed in Ref. [21]. 

The above-mentioned correspondence between the 
Helmholtz and Schrodinger equation is established by 
means of an effective potential [22] that depends not only 
on position and the respective index of refraction n, but 
also on the energy. Changes in the refractive index give 
rise to steps in the effective potential which allows for 
a quantum-mechanical interpretation (e.g., quasibound 
states, tunneling escape). We will discuss this point in 
the next section when we contrast optical systems gov- 
erned by Maxwell's equation with quantum mechanical 
problems obeying the Schrodinger equation. Also, we 
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FIG. 1. (a) Geometry and notations of the dielectric annu- 
lar billiard, (b) Example of a resonant wave function excited 
by a plane wave with wave number k = 6.251 incident from 
the left = 1,R 2 = 0.6,8 = 0.22, m = 3,n 2 = 6). Dark 
regions denote high electric field intensity. 



II. THE DIELECTRIC DISK 

In this section we introduce the methods, techniques, 
and notations used later in the discussion of the annu- 
lar billiard. We present the ray and wave picture for the 
description of optical (or dielectric) systems using the 
simple example of a dielectric disk, which provides all 
the ingredients to deal with the annular billiard (apart 
from a coordinate transformation, see below). We start 
with the ray optics approach and show how methods 
well established in classical dynamics can be adopted 
to optical systems. In the wave description we distin- 
guish between an approach to the resonant states of the 
(naturally) open optical system by complex wave vectors 
based on Maxwell's equations on the one hand, and by 
real wave vectors arising in an 5*-matrix approach on the 
other hand. 



will see how the two possible polarization directions af- 
fect the Maxwell-Schrodinger correspondence and which 
quantity takes the role of h: Maxwell's equations are, of 
course, not aware of the existence of h. 

The further outline of the paper is as follows: In Sec. II 
we introduce the ray and wave optics notion for the sim- 
ple system of the dielectric disk that arises from the annu- 
lar system upon removal of the inner obstacle (ni = n 2 , 
or i?2 = 0). We describe the methods used for the study 
of the annular billiard in the subsequent sections; namely 
the adaptation of the Poincare surface of section method, 
well-known from classical mechanics, to optical systems, 
the exact solution of the Maxwell equation (leading to 
an effective Schrodinger equation) , and the S- matrix ap- 
proach. Whereas the first approach is based on the ray 
picture, the latter two clearly fully include the wave na- 
ture of light. We employ these methods for the annu- 
lar billiard in Sec. Ill, where we introduce methods to 
study the ray dynamics in optical compound systems and 
apply, for the first time to our knowledge, an S'-matrix 
formalism to optical billiards. The expected ray-wave, 
or classical-quantum, correspondence is established in 
Sec. IV and investigated from various viewpoints, includ- 
ing both real space and phase space arguments. However, 
several features in the behaviour of waves require im- 
provements of the simple ray model as we will illustrate 
and explain with typical examples. In our conclusion, 
Sec. V, we discuss the possibility of an experimental re- 
alization of the annular system with the currently avail- 
able dielectric materials. The successes of the ray-picture 
illuminated here and elsewhere [5,6,10,12,23] suggest the 
ray-based design of micro-optical cavities for, e.g., future 
communication technologies. 



A. Ray Optics: Classical billiards with total internal 
reflection 

Within ray optics, the zero-wavelength limit of wave 
optics, light is described by a ray that follows a straight 
line through a medium, very similar to the dynamics of 
a point mass. Let us assume a light ray, or plane wave, 
incident under an angle xi with respect to the normal of 
a dielectric boundary where the refractive index changes 
from ni to n-i- At the interface, the ray is i) specularly re- 
flected under an angle xi = Xi = X: with a polarization- 
dependent [24] probability R™/ TE y see Fig. 2(b). The 
remaining part, T™/ TE = 1 - R™/ TE 7 i s ii) transmit- 
ted into the other medium under an output angle r\ given 
by Snell's law, sin 77 = (ni/712) sinx = nsinx- In the last 
identity we have employed the scaling properties of the 
system that allow one to fix one of the refractive indices, 
(e.g. that of the environment) to unity without loss of 
generality. 

Snell's law cannot be fulfilled to yield real rj for any 
angle of incidence x if n > 1 (ni > 77,2, respectively). 
Total internal reflection occurs if sinx > sinxc = 1/n 
where we introduced the critical angle Xc- For angles of 
incidence above the critical angle, light is confined by to- 
tal internal reflection with zero transmission and behaves 
like a classical point particle, Fig. 2(a). Therefore, real 
and phase space methods from classical mechanics (such 
as ray tracing or the Poincare surface of section tech- 
nique) , prove to be very useful if they are complemented 
by the optical property of refraction: The Poincare sur- 
face of section (SOS) method works exact (except for the 
exponentially small tunneling losses) as long as we are 
in the regime of total internal reflection. However, for 
— 1/n < sinx < 1/n, light can escape so that the in- 
tensity remaining inside the disk is diminishing. This 
fact has to be taken into account when discussing the 
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FIG. 2. Comparison of ray and wave picture for a di- 
electric disk of refractive index n. The upper panels il- 
lustrate the two possibilities of (a) total internal reflection, 
sinx > sin^c = l/n, when the ray dynamics is equal to that 
of a classical point particle in a closed system, and (b) ray re- 
fraction when, due to partial transmission, the light intensity 
inside the disk decreases with time. In (c) the intensity of the 
electric field (see Sec. II B, higher intensity in darker regions) 
for a quasibound state (nkR = 11.428 — 0.254i) of the dielec- 
tric disk (n — 3) is shown. For comparison, an eigenstate of 
the closed disk (nkR — 9.761, vanishing intensity outside the 
disk) is given in the lower left sector (both are for transverse 
magnetic field, TM polarization). 




FIG. 3. Poincare surface of section for the annular billiard 
with Ri = 1,J?2 = 0.6, and S = 0.01. The horizontal axis 
is the polar angle <j>, the vertical axis, sinx, i s proportional 
to the angular momentum in ^-direction (perpendicular to 
the system plane). Although the displacement of the inner 
disk is rather small, it has major impact on trajectories that 
explore the region sinx < R2 + S. Other trajectories are not 
influenced and identical to those of a single dielectric disk, 
see e.g. the two upper-most examples on the right. The 
initial angular momentum of the trajectories shown is positiv 
(sinx > 0). Nonetheless regions where sinx < are explored, 
implying a change in the sense of rotation as, e.g., in the 
lower-most trajectory on the right. The critical angle for total 
internal reflection is marked on the left indicating that the two 
lower-most trajectories are not confined. 



Poincare SOS for optical systems. Figure 3 shows an ex- 
ample of a Poincare SOS for a hard-wall annular system 
with slightly eccentric inner disk (<5 = 0.01). The critical 
value sin x = l/n is marked by an arrow. 

Probing the phase space structure of a rotational in- 
variant system like a disk in terms of a Poincare SOS 
gives a uniform structure as shown in the upper part of 
Fig. 3. Although this is a Poincare SOS for an annu- 
lar system [25] with a slightly eccentric inner disk (see 
Sec. IV), it is identical to that for a disk for trajecto- 
ries that do not hit the inner disk, i.e. sinx > R2 + 8 
(Ri = 1). The straight horizontal lines directly express 
the conservation of angular momentum, that is, conser- 
vation of sinx, and the corresponding trajectories are 
referred to as whispering gallery (WG) orbits. 

The reflection and transmission probabilities, R™/ TE 
and T™/ TE , are provided by Fresnel's laws [26]. A plane 
electromagnetic wave incident on a planar dielectric in- 
terface with angle of incidence x 1S reflected with the 
polarization-dependent probabilities 

r tm = sin 2 (x-??) ^te = tan 2 (x-??) ^ 
sin 2 (x + ?7) ' tan 2 (x + ?y) ' 



where TM (TE) denotes transverse polarization of the 
magnetic (electric) field at the interface, and i] = 
arcsin(n sinx) is the direction of the refracted beam ac- 
cording to Snell's law. 

B. Wave Picture: From Maxwell to Schrodinger 

Before we turn to the more complicated annular bil- 
liard in Sec. Ill, we assume an infinite dielectric cylinder 
of radius R and refractive index n embedded in vacuum 
with refractive index no = 1. We will call k the wave 
number outside the cylinder and, analogously, nk is the 
wave number inside. The solution of Maxwell's equa- 
tions for the vortices of the electromagnetic field [26] is 
given, e.g., in Refs. [4,22], and leads to an equation for 
the electric (magnetic) field that is similar to the conven- 
tional Schrodinger equation. The vector character of the 
fields implies, however, that one has to distinguish two 
possible polarization directions with differing boundary 
conditions. The situation where the electric (magnetic) 
field is parallel to the cylinder (z-) axis is called TM 
(TE) polarization, with the magnetic (electric) field be- 
ing thus transverse. Using the rotational invariance of 
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the system, separation in cylindrical variables (assuming 
a ^-dependence e %m<t ' , and a z-dependence e lkzZ ) eventu- 
ally leads to an effective Schrodinger equation [27] for the 
radial component of the electric field, 



j^TM/TE _ ex p^ n I m ^_R) COSX] . 



(5) 



dr 2 



ld_ 

r dr 



E(r) + V eS {r)E(r) = fc 2 £(r), (2) 



where we introduced the effective potential 

V cS (r) = k 2 (l-n 2 )+m 2 /r 2 + k 2 z . (3) 

The first term reveals immediately that dielectric regions 
with n > 1 correspond to an attractive well in the quan- 
tum analogy, and that a potential structure is determined 
by the change of the refractive indices for different re- 
gions. Note, however, the energy-dependent prefactor - 
a far-reaching difference in comparison to quantum me- 
chanics. The other two terms in Eq. (3) arise from the 
conservation of the angular momentum (characterized by 
the quantum number m), and express the conservation 
of the linear momentum along the cylinder axis (acting 
as an offset in energy), respectively. 

In the following we will consider a dielectric disk (that 
is, we choose a particular cross sectional plane of the 
cylinder to obtain an effective system), and set k z to zero 
corresponding to a wave in the x-y plane. Approaching 
the disk from the outside (r > R, no = 1) there is only 
the angular momentum contribution to the effective po- 
tential V e B, Eq. (3). At r = R, there is a discontinu- 
ity in V e fi that is proportional to 1 — n 2 , reflecting the 
non-continuous change in the refractive index. It reaches 
from fc 2 lax = m 2 /R 2 to k 2 lin = m 2 /{nR) 2 . Inside the 
disk the angular momentum contribution, now shifted 
by fc^ax — fc^im again determines the behaviour (see also 
Fig. 9 at values r/R\ > R2/R1 )• 

The form of the potential suggests an interpretation in 
the spirit of quantum mechanics with metastable states 
in the potential well that decay by tunneling escape, and 
indeed this turns out to be the quantum-mechanical ver- 
sion of confinement by total internal reflection [4,22]. To 
this end we employ a semiclassical quantization condi- 
tion for the z-component of the quantum mechanical and 
classical angular momentum, rah — nhRe(kR) sin \. We 
find 



sinx 



m 



nRe(kR) 



(4) 



as relation between the angle of incidence as ray picture 
quantity, and the wave number and angular momentum 
of a resonance. 

Another correspondence between ray and wave quan- 
tities exists between the (polarization-dependent) Fres- 
nel reflection coefficient Ji™/ TE anc j the imaginary part 
of the wave number that describes the decay of a reso- 
nant state. In fact one can deduce a reflection coefficient 
r tm/te ofthe digk r 28 . 



We wish to point out that there exist deviations between 
the Fresnel values _R™/ TE and ^™/ TE w hen the wave- 
length becomes comparable to the system size, in partic- 
ular around the critical angle. This can be understood 
within a semiclassical picture based on the Goos-Hanchen 
effect [29,30]. 

The general solutions of the radial Schrodinger equa- 
tion (2) are Bessel and Neumann functions, J m (kir) and 
Y m (kir) of order m, where ki is the wave number in the 
respective medium. Since physics requires a finite value 
ofthe wave function at the disk center, the solution inside 
the disk can consist of Bessel functions only. Outside the 
dielectric we assume an outgoing wave function, namely 
a Hankel function of the first kind, in accordance 
with our picture of a decaying state. The resonant states 
are obtained by matching the wave field oc J m {nkr)e lm<t ' 
inside the disk at r = R to the wave field oc H&\kr)e im ' t ' 
outside the disk according to the polarization dependent 
matching conditions deduced from Maxwell's equations. 
The resonant states are solutions of 

J m (nkR) H^'(kR) = V J' m {nkR) ff« (kR) , (6) 

where V = n (1/n) for TM (TE) polarization (primes 
denote derivatives with respect to the full arguments nkr 
and kr, respectively). 

One example of a quasibound state as solution of the 
optically open system is shown in Fig. 2(c) and compared 
to a solution of the closed disk. The shift in the wave pat- 
terns is clearly visible. Owing to the symmetry of the sys- 
tem we find a characteristic (quantum-mechanical) node 
structure that is directly related to the quantum numbers 
m (there are 2m azimuthal nodal points) , and p counting 
the radial nodes (hence m = 3, p = 2 in the example). 

At this point a further discussion concerning the ap- 
pearance of h in optical systems is convenient. Em- 
ploying the quantum-classical correspondence, one ex- 
pects Ti to be related to the reciprocal wave number, 
h ~ 1/fc, because fi — > in the classical (here the ray) 
limit k — > 00. This relation is indeed obtained when 
we compare Eq. (2), divided by k 2 (thereby removing 
the energy-dependence of the effective potential), with 
Schrodinger's equation, and identify 1/k with h. 



C. S-Matrix approach to the dielectric disk 

The main idea when considering a scattering problem 
is to probe the response of the system to incoming (test) 
waves, and to extract system properties like resonance 
positions and widths from the scattered wave. Physically, 
this method is formulated for real wave vectors. 

Here we want to investigate the scattering properties of 
the dielectric disk for electromagnetic waves in the frame- 
work of S'-matrix theory [31-33]. One possible choice for 
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the incident test waves are, of course, plane waves. For 
our rotational invariant disk of finite dimension, however, 
incident waves that allow for angular momentum classifi- 
cation are much more convenient: Then we need to take 
into consideration only waves with impact parameter of 
the order of the system dimension or smaller. The Han- 

(2) 

kcl functions H, n of the second kind possess the desired 
properties. 

Again, we consider a dielectric disk of radius R and 
refractive index n and denote the vacuum wave number 
by k. According to Maxwell's equations and the discus- 
sion in the previous Sec. II B we write the wave function 
^scatt oytgide t na t ig excited by an incident wave of an- 
gular momentum m as 

oo 

*£ att (Ax) =H^(kr)e"^+ £ S^H^ {kr)e a + . 



Here, S m i is the amplitude for an incident wave H$ to 
be scattered into hP'K The scattering amplitudes are 
comprised in the S-matrix. It follows from flux conserva- 
tion that S has to be unitary, a property that we will use 
subsequently. Starting with a general situation in which 
S can have entries everywhere, symmetry requirements 
will reduce the number of independent matrix elements. 
For the dielectric disk, the scattered wave has to obey 
angular momentum conservation and will, therefore, be 
a Hankel function of the same order m as the incoming 
Hankel function. Hence the scattering matrix is diago- 
nal. In the general case where angular momentum is not 
conserved (as, for example, for a deformed disk or the 
annular geometry that we will consider in all following 
sections), scattering will occur into all possible angular 
momenta I. 

Employing the matching conditions (cf. Sec. II B) for 
TM polarization, we obtain from the requirement of con- 
tinuity of the wave function (or the electric field) and 
their derivative eventually the matrix elements S mm ' ■ 



T {k)=A^kT W {k l ) 



(9) 



H^ Y {kr)-n 



JLinkr) rr(2) 



tkr) 



H%>(kr) 



H rn (kr)-n^^rH^(kr 



(7) 



The general idea for identifying resonances is that a 
probing wave with resonance energy will interact longer 
with the system than a wave with "non-fitting" energy 
This can be quantified in terms of the Wigner delay time 
t w (Ek) [34] that is the derivative of the total phase of 
the determinant of the S- matrix, detS=e tS , with respect 
to energy Ek = fc 2 , 



(Ek) 



d6(E k ) 
dEk 



(8) 



In the following, we will use the wave-number based 
delay time 



in order to identify resonances as depicted in Fig. 4(a). 
The solid line shows the result for a dielectric disk with 
n=3. Families of whispering gallery modes (WGMs) are 
identified upon increasing the wave number and can be 
labelled by the quantum number m that counts the az- 
imuthal nodes (2m) . The decrease in peak width, accom- 
panied by an increase in height, that is observed with 
increasing m corresponds to an increase of the angle of 
incidence, Eq. (4), and improved confinement by total 
internal reflection. 

Note the relation between the (total) phase 9(Ek) of 
the S-matrix and the so-called resonance counting func- 
tion: N(E k ) = 6{E k )/2-K, cf. [31]. The idea is that a 
resonance is encountered whenever the phase 6 of det S 
increases by 27r upon increasing the energy E k - 

In the following we will use the function r(k) to de- 
termine the resonances. Isolated resonances appear as 
(Lorentzian) peaks in r(fc), see Fig. 4(a), above a small 
background. Information about the imaginary part of 
the resonance is now encoded in the height and width 
of the Lorentzian resonance peaks [32]. We point out 
that the resolution of very broad and extremely narrow 
resonances might be difficult, because they are either in- 
cluded in the background or not captured using a finite 
numerical grid interval. However, resonances with a wide 
range of widths are easily identified, in particular all res- 
onances that are of interest for microlaser applications 
are found within the ^-matrix approach [35]. 

The area under the curve r(fc) is proportional to the 
number of states with wave numbers smaller than k 
[31,33]. In the case of stepwise potentials such as re- 
alized in ray-splitting billiards, simple Weyl formulas for 
the smooth part of the density of states were derived for 
a number of geometries [18]. The application of these 
results to optical systems where ray splitting is realized 
by refraction and transmission at refractive index bound- 
aries is tempting. However, here we work with an energy 
dependent effective potential, in contrast to the situation 
studied in [18] where only a (stepwise) spatial dependence 
of the potential was assumed. Consequently, a general- 
ization of the formulas derived in [18] would be required 
if one is interested in an analytical expression for the 
smooth part of the density of states, which is, however, 
not the subject of this work. 



III. ANNULAR BILLIARD IN THE RAY AND 
WAVE PICTURE 

In this section we adopt the ray and wave methods 
explained above to the general case of the dielectric an- 
nular billiard. We will denote the three different regions, 
namely the environment (refractive index no), the an- 
nular region (ni), and the inner disk (n 2 ) by the in- 
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FIG. 4. (a) Resonances in the concentric annular billiard 
(Ri = 1, R2 = 0.6, ni =3), corresponding to the first family 
of whispering gallery modes. The annular systems ni = 2 
(dotted) and ni = 4 (dashed) are compared with the ho- 
mogeneous disk 712 = 3. Note the systematic deviation of 
the resonance position to larger (smaller) wave numbers for 
ri2 = 2 (712 = 4) that decreases with increasing angular mo- 
mentum quantum number m since the inner disk becomes 
less important. In the inset, we compare the positions and 
widths of the 4th resonance (m — 4) in the delay time with 
the respective complex wave numbers according to Eq. (10). 
We find excellent agreement with the numerically exact val- 
ues k = 2.0108 - 0.0041i (ti 2 =4), k = 2.0753 - 0.0063i 
(n 2 = 3), and ko = 2.1035 - 0.0075i (n 2 = 2). Note the 
existence of additional resonances for 712 = 4, some of them 
marked by arrows from below and illustrated in the ray and 
wave picture in (b). They are due to the double- well structure 
of the effective potential and referred to as "double WGMs" , 
sec text. 



dices 0, 1, and 2, respectively. The corresponding wave 
numbers are ko,ki, and k 2 . Due to the scaling proper- 
ties we fix tiq = 1 and one of the geometry parameters 
i?i,i?2,<5; we choose R\ = 1. Given a set of parameters 
(n = ^,ni,n 2 ), the same results hold for the scaled set 
(noj^i = nino,7i2 = ^2^0) for wave numbers k = k/nQ, 
if the geometry is not changed. In turn, fixing the di- 
electric constants, the parameter sets (i?i = l,R 2 ,S) 
and (i?i,i?2 = R 2 Ri,S = SR±) are equivalent when 
k^k = k/Ri. 



FIG. 5. Poincare SOS taken at the outer boundary for the 
hard-wall annular billiard (R2 = 0.6,5 = 0.22). Trajecto- 
ries are bounded to the annular region, no optical properties 
of the system are yet included. Typical regular and chaotic 
trajectories are shown both in real and phase space. 



A. Ray optics and refractive billiard 

The rotational invariance of the circular billiard dis- 
cussed in Sec. II can be broken either by deformation (as 
in [12]) or by placing off-centered opaque obstacles inside 
the disk, leading to the hard-wall annular billiard. Start- 
ing from the concentric situation, the system stays close 
to integrable due to the existence of adiabatic invariants 
for not too big eccentricities, see Fig. 3. However, in 
general the phase space of the annular billiard is mixed, 
with regular islands placed in the chaotic sea as shown 
in Fig. 5. 

For optical systems, both the outer and inner bound- 
ary become permeable. Leakage at the outer boundary 
occurs for — 1/n < sinx < 1/n. In the simplest qualita- 
tive picture, starting with the hard-wall system, we will 
assume those rays to leave the cavity (thus simplifying 
Fresnel's laws (1) to a stepwise function). The corre- 
sponding trajectories are assumed to not exist in an op- 
tical cavity. If one is interested in how the intensity of a 
certain trajectory decreases, Fresnel's laws can easily be 
taken into account accurately, leading to the model of a 
Fresnel billiard [5,23,36]. 

However, the description of the refractively opened in- 
ner boundary turns out to be rather complicated. There, 
all rays remain in the billiard, causing a tremendous 
increase of the number of rays upon partial reflection. 
Another crucial difference is that now new trajectories 
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FIG. 6. Poincare SOS taken at the outer boundary for the 
refractive billiard with no = l,ru = 3, rii = 1 and the same 
geometry as in Fig. 5. The hard- wall condition of Fig. 5 is kept 
at the outer boundary, but replaced at the inner one by the 
condition of total internal reflection. If it is not fulfilled the 
ray will penetrate the inner disk, giving rise to a restructured 
phase space and new regular orbits as the one on the lower 
right. 
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FIG. 7. Poincare SOS for the refractive billiard with 
no = l,ni = 3, n.2 = 6 and the same geometry as in Fig. 5. 
Rays in the annulus that hit the inner boundary will pene- 
trate it. Note the existence of whispering gallery modes in the 
inner disk (confined by total internal reflection) not visible in 
this SOS. 



arise, namely those crossing the inner disk. We model 
this situation by introducing the model of the "refractive 
billiard": Whenever total internal reflection is violated 
at the inner boundary, the ray enters the inner disk ac- 
cording to Snell's law with full intensity, such that now 
ray splitting occurs. Otherwise, the ray is specularly re- 
flected and stays in the annulus. This corresponds again 
to a stepwise simplification of Fresnel's laws. Note that 
the hard-wall billiard is in fact a realization of constant 
reflection coefficient Jl™/ TE = i. The real situation 
is found in between the stepwise and constant approxi- 
mations and, depending on the refractive indices chosen, 
results from both limits are needed in order to under- 
stand the resonant modes found in the wave picture, see 
Sec. IV. 

We complete our refractive-billiard model by first as- 
suming specular reflection at the outer boundary, and 
discuss outer-boundary losses subsequently as outlined 
above. Results are shown in Figs. 6 and 7 for the same 
geometry as in Fig. 5 [37], and two different combina- 
tions of refractive indices. In Fig. 6, the annular index 
n\ is highest, allowing for total internal reflection at both 
boundaries. In the limit n\ — > oo we would recover the 
phase space of the hard-wall billiard, Fig. 5. For moder- 
ate Tlx = 3 (no = fi2 — 1) as in Fig. 6 we are, however, 
away from this limit: most of the regular trajectories of 
the hard-wall system are gone and, in turn, new regular 
orbits passing through the inner disk appear. 

The situation changes once more for no < n\ < n,2, 
because then total internal reflection at the inner bound- 
ary is never possible (again, we base our discussion on 
rays entering from the annulus), and all rays hitting the 
inner boundary will enter. Furthermore, they will leave 
the inner disk upon the next reflection according to the 
principle of reversibility of the light path. Note, however, 
that confinement by total internal reflection in the inner 



disk is well possible. From our discussion in Sec. II we 
know that these orbits will and can only be whispering 
gallery modes (WGMs). To anticipate results of the next 
section, those modes do exist and leave their signature 
as very sharp peaks in the delay time. 

In Fig. 7 an example of the phase space is given, show- 
ing yet another structure owing to the change in the re- 
fractive indices. For the regular orbits shown at the right, 
we expect only the upper one to survive the (optical) 
opening of the outer boundary as long as no/ni > 3.2. 
The lower orbit hits the outer boundary perpendicular 
(x = 0) at least in some points, and can therefore only 
be confined by hard walls. 



B. Wave picture: Maxwell's equations and S-matrix 
approach 

Generalizing the wave picture approaches presented in 
Sec. II for the dielectric disk to the annular billiard re- 
quires essentially to consider another, off-centered circu- 
lar boundary at which the matching conditions resulting 
from Maxwell's equations have to be fulfilled as well. An 
eccentric inclusion lowers the rotational symmetry of the 
system to axial reflection invariance about the symmetry 
axis of the system. Consequently, angular momentum is 
not conserved, and the S'-matrix of the compound system 
cannot be diagonal in the general case. 

Maxwell's equations can be solved analytically in the 
concentric case (i?2 > 0), and resonant states with com- 
plex wave number are obtained as zeros of the expression 



iJ m (k2R2)H^'{k R 1 ) x (10) 
H^'ik^H^ihR,) - ffW'(fe 1 E 2 ) J ffi 2 )(fc 1 ii 1 ) 

iJ m {k 2 R2)H^(koR 1 ) x 
HW \k x R 2 )HW \kxRx) - H^'ihR^H^'ihRx) 
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Wave number k 



FIG. 8. Resonance peaks for increasing displacement 8 in 
comparison with the concentric case (dotted) for ni > 1%%. 
Here, whispering gallery modes of the outer disk are shifted 
to the left due to an increase of the effective refractive index. 
However, other modes are affected in a different way, as for 
example the resonances marked by arrows. 



n 2 



j' m (k 2 R 2 )H^'(k Rx) 



x (kiR 2 )H^ (kiRi) - h£> (kxR 2 )H$ (kxRx) 

n in2 j' m {k 2 R 2 )H^ '{koRi) x 

H^{kxR 2 )H^'(kxRx) - Hg\kxR 2 )H^' '(kxRx) 



for TM polarized light. Note that Eq. (10) reduces to 
Eq. (6) for m = n 2 when the annular billiard is reduced 
to a disk. 

In order to investigate the eccentric case, we focus on 
the S-matrix method. The derivation of the S'-matrix 
for the eccentric annular billiard is outlined in the ap- 
pendix. As discussed in Sec. II the information on reso- 
nance position and width is contained both in the com- 
plex wave vector that solves the resonance equation de- 
duced from Maxwell's equations and in the delay-time 
plot r(k). This is illustrated in the inset of Fig. 4(a) 
where the resonance positions and widths found from 
r(fc) are compared with the numerically exact solutions 
of Eq. (10) for concentric geometries. The delay-time plot 
in Fig. 4(a) reveals a systematic deviation of the first few 
resonance positions to the right (left), if the refractive 
index of the inner disk is lower (higher) than that in the 
annulus. However, the deviation from the concentric case 
is rather small. It suggests that the low-lying resonances 
in the (concentric) annular geometry are very similar to 
the WGMs of the dielectric disk and mainly localized at 
the outer boundary. However, the resonant wave func- 
tion does experience the change of the refractive index in 
the inner disk as indicated by the shift of the resonance 
position. The direction of the shift is most easily seen 
when thinking in terms of an effective refractive index 
n G ff, 



An inner disk of lower refractive implies n e g < n\ and a 
larger spacing between the resonances. This is easily un- 
derstood when considering an eigenvalue nk =const. of 
the (closed) dielectric disk. Obtaining the same constant 
value for a smaller n requires a higher k. In contrast, an 
inner disk of higher refractive index reduces the spacing 
between the resonances. This effect is strongest for res- 
onances of high radial quantum number p and small an- 
gular momentum quantum number m, since they rather 
extend to the inner regions of the disk or the annular 
billiard. (In terms of the ray picture, they correspond to 
smaller angles of incidence, leading to the same conclu- 
sion.) Accordingly, the effect reduces for increasing m 
and eventually vanishes if the inner disk is not seen any 
more [38]. In Fig. 4(a) resonances are marked by arrows 
that exist only if the refractive index of the inner disk is 
highest. One corresponding wave pattern, together with 
a ray analogue, is shown in Fig. 4(b). It reveals that the 
"double WGM" structure results from a star-like trajec- 
tory. 

In Fig. 8 we consider the same refractive indices (i.e. 
Tlx = 3,n 2 = 4) and shift now the inner disk off the 
centre. The "double WGMs" (again marked by arrows) 
are affected in a way different from the "conventional" 
WGMs. First of all, the systematic shift of the latter 
can again be understood in terms of the effective refrac- 
tive index. The impact of an off-centered (inner) disk 
is enhanced because in the constricted region it acts like 
a concentric disk with larger radius R^ > R 2 . Note 
that the resonances marked by arrows in Fig. 8 change 
their character from "double WGMs" in the concentric 
and slightly eccentric cases to "generalized" WGMs sim- 
ilar to the one shown in the right panel of Fig. 9 if the 
symmetry breaking caused by the off-centered inner disk 
becomes too strong. 



IV. RAY- WAVE-CORRESPONDENCE FOR THE 
ANNULAR BILLIARD 

In the previous sections we already referred to the ray- 
wave correspondence in optical systems and gave several 
examples which were mainly based on WG modes, and 
the concentric annular billiard. In this section, we first 
continue with WGMs and show how they can be specifi- 
cally influenced by choosing appropriate materials. How- 
ever, ray-wave correspondence holds for far more inter- 
esting trajectories, and we will give illustrative examples 
how closed-billiard trajectories are recovered in the open 
system using real and phase space portraits. 
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FIG. 9. Effective potential for the annular billiard for two 
different sets of refractive indices (no=l). In the lower panels 
corresponding examples of ray trajectories (left) and wave 
functions (right) are shown. The similarity of these reso- 
nances to whispering gallery modes suggest their classification 
as "generalized" WGMs. 



A. Classes of whispering gallery modes in annular 
systems 

In Sec. II we introduced the concept of the effective po- 
tential as a wave-picture method when we established the 
analogy between Helmholtz and Schrodinger equations. 
The generalization of this concept to the annular billiard 
is straightforward and (in the concentric case) essentially 
given by the superposition of two disks, see Eq. (3). The 
result is schematically shown in Fig. 9. Again, we have 
to distinguish two cases: In Fig. 9(a), the refractive in- 
dex rii in the annulus is highest, whereas in Fig. 9(b) the 
optical density increases towards the inner disk. Con- 
sequently, in the first case, Fig. 9(a), the potential well 
coincides with the annular region: Rays in between the 
two disks can be totally reflected at either boundary as 
illustrated in the lower panels. 

The situation is different when the refractive index ni 
of the inner disk is highest [Fig. 9(b)]. The well extends 
now beyond the inner boundary to values r < i? 2 indi- 
cating that the inner disk may support annular WGMs in 
the constricted region, see the lower panels. This is con- 
sistent with the ray picture interpretation stating that 
each ray in the annulus that hits the inner boundary will 
enter the inner disk. However, because of the double-well 
structure of the effective potential this case is even richer: 
There are modes that mainly live in one of the two wells, 
corresponding to WGMs of the inner and outer disk, re- 
spectively. The height of the separating barriers depends 
on the wave number, the quantum number m, the ge- 
ometry, and in particular on the ratio of the refractive 
indices that can be used to tune the height of the barrier 
(note that at the same time the depths of the minima are 
changed) . 



B. Towards closed systems 

Varying the refractive index of an optical system allows 
one to describe the transition between closed and open, 
optical, systems as mentioned earlier. To illustrate this 
fact we increase the refractive index of the annular region 
that we assume to be embedded in vacuum (uq = 712 = 1, 
wave vector fco = A; 2 ). The length spectrum, or Fourier 
transform, of the delay time r(ko) is shown in Fig. 10(a) 
for refractive indices n\ = 3 and n\ = 6 (dashed and full 
line, respectively). The Fourier analysis is performed in 
the spirit of trace formulas that provide a semiclassical 
interpretation of quantum-mechanical results in terms of 
classical periodic orbits for quantum billiards. The quan- 
titative extension of this approach to optical systems will 
require further discussion. Here, we are only interested 
in a qualitative interpretation. 

We have divided the optical length that results from 
the Fourier transformation by n\ in order to compare 
both spectra in terms of geometrical lengths L. The peaks 
in both spectra are rather broad and correspond roughly 
to the circumference of the bigger disk (and higher har- 
monics) which indeed is a typical trajectory length in 
this geometry not only for WGMs, but also for the tra- 
jectory examples shown Figs. 5 and 11 (trajectory parts 
in the inner disk will contribute a length that has to be 
corrected by a factor n 2 /ni). However, the length spec- 
trum for ni — 6 shows an additional peak (marked by the 
arrow) at higher L. A ray trajectory of suitable length 
is shown in Fig. 10(b). The Poincare fingerprint of this 
orbit, see Fig. 5, possesses regular islands at sinx = 0, 
where in the simplest interpretation refractive escape will 
occur, independent on the refractive index! That modes 
of this type are found for sufficiently large n\, indicates 
that we have to refine our interpretation. For example, 
we can discuss the Fresnel reflection coefficient R± at 
normal incidence, R± — {tlq — ni) 2 /(no + ni) 2 that in- 
creases as Ui is increased, reaching the value one in the 
limit ni — > 00, in accordance with the picture of complete 
internal reflection. This explains the observed behaviour 
and we present more examples of "sophisticated" ray- 
wave, or classical-quantum, correspondence in the next 
paragraph. 

C. Correspondence in real and phase space 

In Fig. 11(a) we show a typical delay-time plot r(fco) 
for the annular billiard with the same geometry as before 
and refractive indices ri\ = 3, 77,2 = 6. We investigate 
low-lying resonances that show a characteristic grouping 
of four resonances (marked by A, B, C, D) over several 
periods. The corresponding wave patterns, together with 
suggestions for ray analogues, are shown in Fig. 11(b) for 
each resonance. We have mainly chosen regular orbits as 
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FIG. 10. (a) Fourier transform of the delay time r(fco) for 
the annular billiard geometry R2 = 0.6, S — 0.22 and two re- 
fractive index combinations, no = U2 = 1, rii = 3 (dashed 
line), and ri\ — 6 (full line). The appearance of a new peak at 
larger geometric length is clearly visible. A suitable (quasi-) 
periodic orbit candidate together with a resonant state is 
shown in (b). 
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FIG. 11. (a) Delay time r(fco) for the annular billiard ge- 
ometry as in Figs. 5 and 10 and ni = 3, 712 = 6. The structure 
is dominated by groups of four resonances labelled A, B, C, 
D. The resonance at fcrj ~ 6.251 is the one shown in Fig. 1(b). 
The four resonances of the second group (6.6 < ftrj < 6.9) are 
shown in (b) together with a ray trajectory representative. In 
part (c), rays (crosses) and waves (intensity plot, high inten- 
sity in dark regions) are compared in phase space in terms of 
their Poincare and Husimi representation, respectively. 



candidates because of the regular structure of the delay- 
time plot. Neighbouring resonances of the same kind 



(i.e., the same letter) indeed differ by one in the number 
of nodes [39] . Note that the ray representatives stem from 
both the hard-wall and the refractive billiard simulation, 
see Figs. 5 and 7. 

In Fig. 11(c) we computed the Husimi function 
[16,20,40] for each of the wave functions [41], and marked 
the rays by crosses in the corresponding Poincare SOS 
such that we can directly compare the phase space pre- 
sentations of waves and rays. The coincidence between 
regular islands and high-probability regions (dark) of the 
Husimi function appears satisfying on first sight. How- 
ever, closer inspection reveals differences in the details. 
For example, Husimi "islands" are shifted away from reg- 
ular islands as in the case of resonance D, with the cor- 
responding real space modifications [Fig. 11(b)] clearly 
visible as well. One possible explanation might be pro- 
vided by the Goos-Hanchen effect that causes a lateral 
shift of the reflected ray for angles of incidence around 
and greater than the critical angle [29,30], thereby effec- 
tively changing the angle of incidence. Furthermore, we 
point out that the ray trajectory for resonance D is known 
from the hard-wall system. The qualitative similarity to 
the corresponding resonant wave pattern is remarkable, 
and one might think of the differences as in order to meet 
the new interference requirements caused by the optical 
opening of the inner disk. This gives yet another exam- 
ple of the predictive power of the simple ray model when 
only the qualitative character of the resonances is of inter- 
est and importance. On the other hand, it proves to be 
essential to consult wave methods when one is interested 
in details. 

V. CONCLUSIONS 

To conclude, we have investigated the ray and wave 
properties of composite optical systems by applying 
methods known from the classical and quantum theory 
of mixed dynamical systems. Using the optical annu- 
lar billiard as an example, we have shown this concept 
to be very fruitful. This means in particular that al- 
ready the simple ray model provides a good qualitative 
understanding of the system properties, even for small 
wave numbers below nkR « 30. However, care must 
be taken when quantitative results are required, or the 
classical (ray) phase space is directly translated into ex- 
pected wave patterns: We find regular orbits associated 
with regular islands in phase space to be the dominant 
class of resonant wave patterns, and suppression of wave 
functions hosted by the chaotic part of the phase space. 
The dependence of this behaviour on the size of the wave 
number (i.e. 1/K) remains an interesting topic for future 
work. 

One remark is due concerning the refractive indices em- 
ployed in the calculations. The index n — 3 often used 
here is higher than that of water (1.33) or glass (around 



10 



1.5 up to 1.8) but is easily reached in semiconductor com- 
pounds where typically n = 3.3. An index n = 6 seems 
to be presently out of reach, which, however, does not 
affect the conclusions drawn here. 

Summerizing, ray picture results may serve as a guide 
in the investigation of wave properties of optical systems, 
even away from the ray limit k — > oo. For the annular 
billiard as an example of a compound cavity system we 
demonstrated that the dominant resonant wave patterns 
can be seen as originating from the regular orbits of both 
the hard-wall and the refractive billiard. This knowl- 
edge can be more generally used, e.g., in the construc- 
tion of microlasers with designed properties. Knowing 
the potential reflection points and high-intensity regions 
of modes from simple ray-based considerations allows one 
to design microcavities with custemized properties. Pre- 
dictions can be made concerning, e.g., the effective cou- 
pling between and into cavities or how to efficiently pump 
lasing systems. In turn, one can think of cavity shapes 
designed according to the technical requirements. The 
application of the ray-wave correspondence in sophisti- 
cated optical (compound) systems therefore may provide 
a powerful tool for future optical communication tech- 
nologies. 

Optical cavities represent interesting model systems for 
quantum-chaos motivated studies. We successfully ap- 
plied the S-matrix approach to gain spectral information, 
and qualitatively discussed its periodic-orbit interpreta- 
tion [43]. The development of quantitative semiclassical 
theories in the spirit of the Weyl and the trace formu- 
las remains an open subject, in particular for compound 
systems consisting of more than one region with fixed 
refractive index like the annular billiard. 



ACKNOWLEDGMENTS 

We thank J. U. Nockel for an introduction into the 
subject of optical cavities and acknowledge many use- 
ful discussions with T. Dittrich, S. Fishman, G. Hacken- 
broich, J. U. Nockel, H. Schomerus, H. Schanz, R. Schu- 
bert, P. Schlagheck, U. Smilansky, and J. Wiersig. M. H. 
thanks U. Smilansky for his hospitality at the Weizmann 
Institute. 




FIG. 12. Addition theorem for Bessel functions. 



Sec. IIC, the situation we are confronted with here is 
more complicated: the two disks lay one in the other, 
and their centres will in general not coincide. 

We will begin with the scattering problem at the inner 
boundary and express the S'-matrix S l of the dielectric 
disk with respect to a coordinate system with origin dis- 
placed from the centre of the inner disk. This implies 
that S % is not diagonal. From Sec. II C we already know 
the (diagonal) S'-matrix S lc of the inner disk in primed 
coordinates (see Fig. 1), Eq. (7). We will now derive the 
relation between S lc and S l . 

To this end we write the ansatz for the wave function 
in the annulus in primed coordinates, r' = r — 5, with 5 
being the vector from the centre of the large disk to the 
centre of the smaller disk, as 



* lc (r-d)= £ af \Hf\h\r-5\) 



I — — OO 



J2 SftHMfair-fye*' 



, (Al) 



where the coefficients af are to be chosen to yield the 
desired kind of incident wave. We use the addition theo- 
rem for Bessel functions Z m G { J m , Y m , Hm\ ffm'} [42] 
to relate the arguments wr to wr' (w is a constant factor, 
and we assumei?2 > $), see Fig. 12, 



Z m {wr')e lm ' 1 ' = ^ Jk(w6)Z m+k (wr)t 



i(m-\-k)(j) 



(A2) 



k=—oo 



APPENDIX A: S'-MATRIX FOR THE ANNULAR 
BILLIARD 

We will generalize the ideas developed in Sec. II C to 
the dielectric annular billiard in order to determine the 
S-matrix for the eccentric annular billiard. This problem 
can be divided into the scattering problem at the outer 
boundary (between refractive indices no and n±) and that 
at the inner boundary (between indices n\ and 712). Al- 
though the scattering at a dielectric disk was solved in 



Inserting this into Eq. (Al), we obtain the expression 



tt 1 



E E 

(— — 00 k=— oc 



xJ fe (M) e i{l+k ^ 



(A3) 



for the wave function in the annulus, now expressed with 
respect to the centre of the larger disk, i.e., in unprimed 
coordinates. We specify the coefficients af by the re- 
quirement that the amplitude of an incident wave with 
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angular momentum m shall be normalized to one in un- 
primed coordinates, 



oo oo 



I — — 00 k— — oo 

With n = l + k, and J2k Jm-{n-k)Jk = <W we find that 
choosing 



l af = J m -l(kiS) VI 



(A4) 



provides a suitable set of coefficients for a given m. Ac- 
cordingly, we write 

oo ( oo 

E E [^Hf\ kl r)e^ 

m= — oo I, fi,k=—oc 

+J m -( AJ - fe )(fci<5)^_ fc)(M _ fc) J fc (fc 1( 5) J ff( 1 )(fc 1 r)e^ 

OO f OO 

m— — oo ^ /j.— — oc J 

where we have read off the scattering matrix S l of the 
inner disk with respect to the centre of the outer disk, 

oo 

Si n , = E J m -(»-k)(ki5)Sfr_ k){tl _ k) J k (hS) . (A5) 

k— oo 

The structure of this equation suggests a notation in 
terms of a transformation matrix U, namely S l = 
U~ l S lc U , that describes the change in the origin of the 
coordinate system. We find Uyi = Jp-i and U^ 1 = Ji-v ■ 
The scattering matrix S l allows us to describe the scat- 
tering at an off-centred disk, and we can now formulate 
the scattering problem of the annular billiard in the spirit 
of Sec. II C. Accordingly, we start with an ansatz for the 
wave function outside the annular system (\r\ > 
using polar coordinates) of the form 

oo 

*V) = E 



M=— oo 



E 

M= — oo 



M'=-oo 



where we have introduced the scattering matrix S of the 
(compound) system and the definitions 



y M (k r) = H$(k ry M * 
*+(fc o r) = ^ ) (fc o r)e 4M0 



(A6) 
(A7) 



for incoming and outgoing waves outside the disk. Note 
that we have used the freedom in fixing one of the am- 
plitudes. 

Similarly, we write for the wave function ^P 1 in the 
annular region 



at 



tfj-(fcif) 



OC- 

E 



S?,,*+(fcirO 



(A8) 



with the amplitudes ai, the abbreviations as in Eqs. (A6, 
A7), and S l from Eq. (A5). 

Now, we determine S from the matching conditions, 
introduce the notation of capital letters for functions of 
argument kor, and reserve lower case characters for the 
argument k\r. Given an incident wave of angular mo- 
mentum M, wave function matching for each angular 
momentum L of the scattered waves yields 



a <M) h W e iL* + J- al M) Si L h¥e iL «, 



AM) 



are coefficients associated 



where the amplitudes Oj 
with an incoming function of angular momentum M, 
namely H M '. Since this has to hold for all M, and at 
fixed M for all L, we write this as a matrix equation 

( (M) jff (2)| + (5 (M)| jff (l) = (a (M) | ^(2) + ^^(l)^ ( (AQ) 

where S l is a matrix, h,W and /i« are diagonal matri- 



ces, h 



1,2) 



(1 2) 

h\ ' Sij, and we adopt the 6ra-notation for 
quantities that, at fixed M, are transposed vectors and 
gain matrix character once M is varied. With this no- 
tation we immediately write the matching condition for 
the derivatives as 

fc (( (M) tf (2) '| + <S (M) |tf (1) ') 

= (a^\k 1 {h^' + S i h^') . 

From Eq. (A9) we find after substituting F = + 
S*h^ that 



,( M )| 



Introducing furthermore F' — h'^ + S l h' (1> and W — 
F~ 1 F', we write the 5-matrix solution of the problem as 

S = (fcr W-k H' {2) ) (k H' {1) h w) ' . 

This last equation allows us to apply the Wigner-delay- 
time approach to resonances, cf. Sec. II C, and we used 
this method to study resonances of the optical annular 
billiard, cf. Sees. Ill and IV. 

We complete the discussion here with some comments 
on the wave functions. We have not yet given the wave 
function in the inner disk. The ansatz is a sum over 
Bessel functions, 

oo 

* 2 V) = E 6f Ji(fer')e a *' , 

I — — oo 
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where we adopted primed coordinates for convenience. 
The coefficients bi are found from matching with the wave 
function in the annulus at the inner boundary. To this 
end we have to rewrite the annular wave function (A8) 
in terms of primed coordinates by applying the addition 
theorem (A2) for Bessel functions. After straightforward 
algebra we find 



* lc (r') 



E 



St?H, 



(i) 



far') 



where the coefficients a c { are related to the a; by a\ = 

Another remark is in order concerning the validity of 
the addition-theorem-based expansion of the Bessel func- 
tion when changing between primed and unprimcd co- 
ordinates. Expansion of the annular wave function in 
primed coordinates fails near the outer boundary where 
\r'\ > i?i — i?2 — <5. Similarly, expanding the annular 
wave function in unprimed coordinates does not work 
near the inner boundary where \r\ < R 2 + S. The rea- 
son for this behaviour is that angular momentum is not 
conserved in the eccentric annular billiard, and the ex- 
pansion breaks down at radii where waves explore this 
symmetry-breaking region because the corresponding in- 
terface boundaries are hit. 
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